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Abstract 

Background: The previous published data on the association between CYP1A2*F (rs762551), CYP1B1 Leu432Val 
(rsl 056836), Asn453Ser (rsl 80040), and Arg48Gly (rsl0012) polymorphisnns and colorectal cancer risk remained 
controversial. 

Methodology/Principal Findings: The purpose of this study is to evaluate the role of CYP1A2*F, CYP1B1 Leu432Val, 
Asn453Ser, and Arg48Gly genotypes in colorectal cancer susceptibility. We performed a meta-analysis on all the eligible 
studies that provided 5,817 cases and 6,544 controls for CYP1A2*F (from 13 studies), 9219 cases and 10406 controls for 
CYPIBI Leu432Val (from 12 studies), 6840 cases and 7761 controls for CYPIBI Asn453Ser (from 8 studies), and 4302 cases 
and 4791 controls for CYPIBI Arg48Gly (from 6 studies). Overall, no significant association was found between CYP1A2*F, 
CYPIBI Leu432Val, Asn453Ser, and Arg48Gly and colorectal cancer risk when all the eligible studies were pooled into the 
meta-analysis. And in the subgroup by ethnicity and source of controls, no evidence of significant association was observed 
in any subgroup analysis. 

Conclusions/Significance: In summary, this meta-analysis indicates that CYP1A2*F, CYPIBI Leu432Val, Asn453Ser, and 
Arg48Gly polymorphisms do not support an association with colorectal cancer, and further studies are needed to 
investigate the association. In addition, our work also points out the importance of new studies for CYP1 A2*F polymorphism 
in Asians, because high heterogeneity was found (dominant model: /^ = 81.3%; heterozygote model: F = 79.0). 
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Introduction 

Sporadic colorectal cancer (CRC) is considered to be a 
multifactorial disease, in which multiple exposures to endogenous 
factors and dietary carcinogens interact with individual genetic 
background in a complex manner resulting in modulation of the 
risk [1]. In 2010, an estimated 142,570 new cases will be diagnosed 
and 51,370 deaths wiU occur in the whole world [2]. Epidemi- 
ologic studies on Western populations have emphasized the large 
contribution of food and lifestyle to sporadic CRC risk [3-7]. 
High-fat and low-fiber diets, as well as alcohol, tobacco, and red or 
processed meat consumption, have been shown to produce high 
levels of polycycUc aromatic hydrocarbons and heterocyclic 
aromatic amines. These procarcinogenic agents are potentially 
very harmful and may play a key role in the malignant 
transformation of cells by interacting with DNA [8] . It has been 
proposed that this risk may be due to carcinogenic polycyclic 



aromatic hydrocarbons (PAHs) and heterocyclic amines produced 
when meat is cooked at high temperatures [9] . 

CYPIBI gene is located on chr2p22-p21, which is involved in 
the metabolic activation of polycyclic aromatic hydrocarbons 
(PAHs) including benzo(a)pyrene and dimethylbenz(a)anthracene 
(DMBA), but with a product distribution that is distinct from 
CYPlAl [10,11]. Several lines of evidence suggest that CYPIBI 
plays a role in carcinogenesis. CYPIBI is commonly over- 
expressed inhumanmalignancies [12] and activates a variety of 
carcinogens. For example, CYPIBI catalyzes both the formation 
of dihydrodiols of specific PAHs and their subsequent oxidation to 
carcinogenic dihydrodiol epoxides [13]. In humans, CYPIBI is 
genetically polymorphic and more than 50 single nucleotide 
polymorphisms (SNPs) have been reported so far, of which certain 
deleterious mutations are associated with primary congenital 
glaucoma [14]. Of the most common SNPs of CYPIBI gene, four 



PLOS ONE I www.plosone.org 



1 



August 2014 | Volunne 9 | Issue 8 | el 00487 



Gene Polymorphism and Colorectal Cancer Risk 



have been reported to result in amino acid substitutions including 
Arg by Gly at codon 48 (rsl0012), Leu by Val at codon 432 
(rsl056836) and Asn by Ser at codon 453 (rsl800440). GYP 1A2 is 
an important gene in catalyzing 2- and 4-hydroxylations of 
estrogens [40-42] and metabolism of carcinogens [43—45]. 
GYP1A2*1G, located in the 5'-non-coding promoter region of 
GYP1A2, was reported to be associated with decreased enzyme 
inducibility in Japanese smokers but seems to be very rare [46] . 

To date, a number of molecular epidemiological studies have 
been done to evaluate the association between CYP1A2*F, 
GYPIBI Leu432Val, Asn453Ser, and Arg48Gly polymorphisms 
and colorectal cancer risk in diverse populations [15-29,31,32,34— 
39]. However, the results were inconsistent or even contradictory. 
Therefore, we performed a comprehensive meta-analysis by 
including the most recent and relevant articles to identify statistical 
evidence of the association between CYP1A2*F, GYPIBI 
Leu432Val, Asn453Ser, and Arg48Gly polymorphisms and risk 
of colorectal cancer that have been investigated. Meta-analysis is a 
powerful tool for summarizing the difiFerent studies. It can not only 
overcome the problem of small size and inadequate statistical 
power of genetic studies of complex traits, but also provide more 
reliable results than a single case-control study. 

Materials and Methods 

Identification and eligibility of relevant studies 

A comprehensive literature search was performed using the 
PubMed, GNKI, and Medline database for relevant articles 
published (the last search update was Sep 10, 2013) with the 
following key words "CYP1A2", "CYPIBI", "polymorphism", 
"Variant", or "Mutation", and "Golorectal". In addition, studies 
were identified by a manual search of the reference lists of reviews 
and retrieved studies. We included all the case-control studies and 
cohort studies that investigated the association between GY- 
P1A2*F, CYPIBI Leu432Val, Asn453Ser, and Arg48Gly poly- 
morphisms and colorectal cancer risk with genotyping data. AU 
eligible studies were retrieved, and their bibliographies were 
checked for other relevant publications. 

Inclusion criteria 

The included studies have to meet the following criteria: (1) only 
the case-control studies or cohort studies were considered; (2) 
evaluated the GYP1A2*F, GYPIBI Leu432Val, Asn453Ser, and 
Arg48Gly polymorphisms and the risk of colorectal cancer; (3) the 
genotype distribution of the polymorphism in cases and ctmtrols 
were described in details and the results were expressed as odds 
ratio (OR) and corresponding 95% confidence interval [95"/o GI). 
Major reasons for exclusion of studies were as follows: (1) not for 
cancer research; (2) only case population; (3) duplicate of previous 
publication (When the same patient population was used in several 
publications, only the most recent, largest or complete study was 
included following careful examination). 

Data extraction 

Information was carefully extracted from all eligible studies 
independentiy by two investigators according to the inclusion 
criteria listed above. The following data were collected from each 
study: first author's name, year of publication, country of origin, 
ethnicity, source of controls (populatioii-bascd controls, hospital- 
based controls, and family-based controls), and numbers of cases 
and controls in the GYP1A2*F, GYPIBI Leu432Val, Asn453Ser, 
and Arg48Gly genotypes whenever possible. Ethnicity was 
categorized as "Gaucasian" and "Asian". When one study did 
not state which ethnic groups was included or if it was impossible 



to separate participants according to phenotype, the sample was 
termed as "mixed population". We did not define any minimum 
number of patients to include in this meta-analysis. Articles that 
reported different ethnic groups and different countries or 
locations, we considered them different study samples for each 
category cited above. 

Statistical analysis 

Grude odds ratios (ORs) together with their corresponding 95% 
confidence intervals {95% GIs) were used to assess the strength of 
association between the GYP1A2*F, GYPIBI Leu432Val, 
Asn453Ser, and Arg48Gly polymorphisms and colorectal cancer 
risk. The pooled ORs were performed for dominant model 
(GYP1A2*F: GY + YY vs. GG; GYPIBI Leu432Val: Leu/Val + 
Val/Val vs. Leu/Leu; GYPIBI Asn453Ser: Asn/Ser + Ser/Ser vs. 
Asn/ Asn; GYPIBI Arg48Gly: Arg/ Gly + Gly/ Gly vs. Arg/ Arg), 
recessive model (GYP1A2*F: YY vs. GG + CY; GYPIBI 
Leu432Val: Val/Val vs. Leu/Leu -I- Leu/Val; CYPIBI Asn453- 
Ser: Ser/Ser vs. Asn/Asn + Asn/Ser; CYPIBI Arg48Gly: Gly/ 
Gly vs. Arg/ Arg + Arg/Gly), co-dominant model (GYPl A2*F: YY 
vs. GG and GY vs. GG; GYPIBI Leu432Val: Val/Val vs. Leu/ 
Leu and Leu/Val vs. Leu/Leu; GYPIBI Asn453Ser: Ser/Ser vs. 
Asn/Asn and Asn/Ser vs. Asn/Asn; GYPIBI Arg48Gly: Gly/ Gly 
vs. Arg/ Arg and Arg/Gly vs. Arg/ Arg), and additive model 
(GYP1A2*F: Y vs. C; CYPIBI Asn453Ser: Ser/ Asn; GYPIBI 
Asn453Ser: Ser vs. Asn; GYPIBI Arg48Gly: Gly vs. Arg), 
respectively. Between-study heterogeneity was assessed by calcu- 
lating (5-statistic (Heterogeneity was considered statistically signif- 
icant if P<0.10) [47] and quantified using the value, Venice 
criteria [48] for the test included: "/ <25% represents no 
heterogeneity, = 25-50% represents moderate heterogeneity, 
7^ = 50-75% represents large heterogeneity, and /^>75% repre- 
sents <;xtrcmc lictc'rogx'neity". If results wcyc not lictcrogeiieous, 
the pooled ORs were calculated by the fixed-effec:t model (we used 
the Q-statistic, which represents the magnitude of heterogeneity 
between-studies) [49] . Otherwise, a random effect model was used 
(when the heterogeneity between-studies were significant) [50] . We 
also performed subgroup analysis by ethnicity and source of 
controls were conducted. Moreover, sensitivity analysis was 
performed by excluding a single study each time. We also ranked 
studies according to sample size, and then repeated this meta- 
analysis. Sample size was classified according to a minimum of 200 
participants and those with fewer than 200 participants. The cite 
criteria were previously described [51]. HWE was calculated by 
using the goodness-of-fit test, and deviation was considered when 
P<0.05. Begg's funnel plots [52] and Egger's linear regression test 
[53] were used to assess publication bias. We opted for using 
ethnicity, source of controls, menopausal status, and sample size as 
possible different sources of heterogeneity. AU of the calculations 
were performed using STATA version 10.0 (STATA Corporation, 
College Station, TX). 

Results 

Literature search and meta-analysis databases 

Relevant publications were retrieved and preliminarily 
screened. As shown in Fig. 1, 43 publications were identified, 
among which 6 irrelevant papers were excluded. Thus, 37 
publications were eligible. Among these publications, 14 articles 
were excluded because they were review articles, case reports, and 
other polymorphisms of GYPl A2 and CYPIBI. As summarized in 
Table 1, 23 articles with 39 studies were selected in this meta- 
analysis, including 5,817 cases and 6,544 controls for GYP1A2*F 
(from 13 studies), 9,219 cases and 10,406 controls for CYPIBI 
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Leu432Val (from 12 studies), 6,840 cases and 7,761 controls for 
CYPIBI Asn453Ser (from 8 studies), and 4,302 cases and 4,791 
controls for CYPIBI Arg48Gly (from 6 studies). Among these 
studies, eight were Caucasians, four were Asians, and 1 mixed 
populations for CYP1A2*F. All studies were Caucasians except for 
one study was mixed population for CYPIBI polymorphisms. The 
distribution of genotypes in the controls was consistent with 
Hardy-Weinberg equilibrium in all studies. All of the cases were 
pathologically confirmed. 

Meta-analysis results 

Table 2 lists the main results of the meta-analysis of 
CYP1A2*F polymorphism and colorectal cancer risk. Overall, 
no significant association was found between CYP1A2*F poly- 
morphism and colorectal cancer risk (dominant model: OR — 1 .05, 
95% CI = 0.94-1.18, Ph = 0.010, 7^ = 54.1%; recessive model: 
OR= 1.01, 95% CI = 0.90-1.13, /'h = 0.426, 7^ = 2.0%; homozy- 
gote model: OR = 1.04, 95% CI = 0.93-1.17, Ph = 0.144, 
7^ = 30.0%; heterozygote model: OR = 1.05, 95% CI = 0.94- 
1.17, Ph = 0.023, 7^ = 49.2%; additive model: OR = 1.03, 95% 
CI = 0.9,5-1.11, Ph = 0.026, 7^ = 48.2%, Fig. 2). Signifirant be- 
tween-study heterogeneity was detected. Hence, we performed the 
stratified analyses according to ethnicity and source of controls. In 
the stratified analysis by ethnicity, no significant association was 
found among Caucasians (dominant model: OR= 1.02, 95% 
CI = 0.9,5-1.10, Ph = 0.233, 7^ = 24.6%; reces.sive model: 
OR= 1.06, 95% CI = 0.94-1.20, Ph = 0.387, 7^ = 5.6%; homozy- 
gote model: OR = 1.07, 95% 01 = 0.94-1.21, 7'h = 0.224, 
7^ = 25.6%; heterozygote model: OR=1.01, 95% CI = 0.94- 
1.09, Ph = 0.403, 7^ = 3.5%; additive model: OR= 1.03, 95% 
CI = 0.97-1.08, Ph = 0.157, 7^ = 34.0%, Fig. 3) and Asians 
(recessive model: OR = 0.78, 95% CI = 0.57-1.05, Ph = 0.681, 
7^ = 0.0%; homozygote model: OR=0.91, 95% 01 = 0.49-1.68, 
Ph = 0.076, 7^ = 56.5%; additive model: OR = 0.98, 95% 
CI = 0.69-1.42, Ph = 0.009, f = 74.3%, Fig. 4). In addition, high 



heterogeneity was found among Asians (dominant model: 
7^ = 81.3%; heterozygote model: 7^ = 79.0). When grouped by 
source of control, there was still no evidence of significant 

association. 

Table 2 also lists the main results of the meta-analysis of 
CYPIBI Leu432Val polymorphism and colorectal cancer risk. 
Overall, no significant association was found between CYPIBI 
Leu432Val polymorphism and colorectal cancer susceptibility 
(dominant model: OR= 1.00, 95% CI = 0.94-1.06, Ph = 0.770, 
7^ = 0.0%; recessive model: OR = 1.05, 95% CI = 0.98-1.13, 
Ph = 0.251, 7^ = 20.3%; homozygote model: OR= 1.04, 95% 
CI = 0.96-1.13, Ph = 0.383, 7^ = 6.3%; heterozygote modeh 
OR=0.98, 95% CI = 0.91-1.04, Ph = 0.687, 7^ = 0.0%; additive 
model: OR= 1.02, 95% CI = 0.98-1.06, Ph = 0.498, 7^ = 0.0%). 

Table 2 also lists the main results of the meta-analysis of 
CYPIBI Asn453Ser polymorphism and colorectal cancer risk. 
Overall, no significant association was found between CYPIBI 
Asn453Ser polymorphism and colorectal cancer susceptibility 
(dominant model: OR = 0.97, 95% CI = 0.87-1.08, Ph = 0.053, 
7^ = 49.6%; recessive model: OR = 0.92, 95% CI = 0.76-1.1 1, 
Ph = 0.617, 7^ = 0.0%; homozygote model: OR = 0.92, 95% 
CI = 0.76-1.11, Ph = 0.685, 7^ = 0.0%; heterozygote model: 
OR = 0.97, 95%. CI = 0.86-1. 11, Ph = 0.016, 7^ = 61.8%; additive 
model: OR = 0.97, 95% CI = 0.91-1.03, Ph = 0.135, 7^ = 38.6%). 
Significant between-study heterogeneity was detected. Hence, we 
performed the stratified analysis according to source of controls. 
And in the subgroup analysis by source of controls, there was stiU 
no significant association detected in any genetic model. 

Table 2 also lists the main results of the meta-analysis of 
CYPIBI Arg48Gly polymorphism and colorectal cancer risk. 
Overall, no significant association was found between CYPIBI 
Arg48Gly polymorphism and colorectal cancer susceptibility 
(dominant model: OR = 0.99, 95% CI = 0.91-1.08, Ph = 0.780, 
7^ = 0.0%; recessive model: OR = 1.00, 95% CI = 0.86-1.16, 
Ph = 0.138, 7^ = 40.1%; homozygote model: OR= 1.00, 95% 
CI = 0.86-1.16, Ph = 0.124, 7^ = 42.1%; heterozygote modeh 
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Figure 1. Study flow chart explaining the selection of the 23 eligible articles included in the meta-analysis. 
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Figure 2. Forest plot of CYP1A2*F polymorphism and colorectal cancer risk among overall analysis (additive model). 

doi:10.1371/journal.pone.0100487.g002 



OR=0.99, 95% CI = 0.91-1.08, Ph = 0.989, r = 0.0%; additive 
model: OR = 0.97, 95% CI = 0.91-1.03, Ph = 0.135, 7^ = 38.6%). 

Test of heterogeneity and sensitivity 

There was significant heterogeneity among these studies for 
dominant model comparison (Pj, = 0.008 for CYP1A2*F and 
Pii — 0.053 for CYPIBI Asn453Ser), heterozygote model compar- 
ison (Ph = 0.020 for CYP1A2*F and Ph = 0.016 for CYPIBI 
Asn453Ser) and additive model comparison (Pi-j = 0.022 for 



CYP1A2*F). Then, we assessed the source of heterogeneity by 
ethnicity and source of controls. We found that ethnicity and 
source of controls [data not shown) did not contribute to substantial 
heterogeneity. Sensitivity analysis was conducted to determine 
whether modification of the inclusion criteria of this meta-analysis 
affected the results. Although the sample size for cases and controls 
in all eligible studies ranged from 175 to 2,455, the corresponding 
pooled ORs were not qualitatively altered with or without the 
study of small sample. In addition, a single study involved in the 
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Figure 3. Forest plot of CYP1A2*F polymorphism and colorectal cancer risk among Caucasians (additive model). 

doi:1 0.1 371 /journal.pone.01 00487.g003 



PLOS ONE I www.plosone.org 



6 



August 2014 | Volume 9 | Issue 8 | el 00487 



Gene Polymorphism and Colorectal Cancer Risk 



Study 

Kobayashi [19] 2009 
Yoshida [22] 2007 
Bae [25] 2006 
Chen [26] 2006 
Overall 




.376390 




Odds ratio 



Odds ratio 
(95% CI) 



% Weight 





0.74 ( 0.52,1.05) 25.8 



0.83 ( 0.53.1.30) 22.3 



1.77 ( 1.18,2.66) 24.0 



0.88 ( 0.66,1.18) 27.9 



0.98 ( 0.68,1.42) 100.0 



2.65681 



Figure 4. Forest plot of CYP1A2*F polymorphism and colorectal cancer risk among Asians (additive model). 

doi:10.1371/journal.pone.0100487.g004 



meta-analysis was deleted each time to reflect the influence of 
individual data set to the pooled ORs. The results were also not 
qualitatively altered. 

Publication bias 

Both Begg's fijnnel plot and Egger's test were performed to 
assess the pubUcation bias of literatures. The Egger's test results 
and Begg's funnel plot (Fig. 5, 6) suggested no evidence of 
publication bias in the meta-analysis of CYP1A2T (P = 0.160 for 
dominant model, P = 0.714 for recessive model, P = 0.862 for 
homozygote model; P = 0.248 for heterozygote model; P = 0.462 
for additive model) and Leu432Val (P = 0.749 for dominant 
model, P = 0.864 for recessive model, P = 0.991 for homozygote 
model; P = 0.721 for heterozygote model; P = 0.689 for additive 
model), although possible publication bias was suggested for 
Asn453Ser polymorphism with colorectal cancer risk in additive 
model and recessive model and for Arg48Gly with colorectal 
cancer risk in any genetic model. This might be a limitation for 
meta-analysis of Arg48Gly and Asn453Ser polymorphisms, 
especially those with small sample size, are less likely to be 
published. Figure 7, 8 lists the Duval and Tweedie nonparamet- 
ric "trim and fill" methods funnel plot in additive model and 
recessive model. Adjusting for possible publication bias using the 
Duval and Tweedie nonparametric "trim and fiU" method for 
overall studies, the results did not change between Arg48Gly and 
Asn453Ser polymorphism with colorectal cancer risk. 

Discussion 

CYPIBI is commonly over-expressed inhumanmalignancies 
and activates a variety of carcinogens. For example, CYPIBI 
catalyzes both the formation of dihydrodiols of specific PAHs and 
their subsequent oxidation to carcinogenic dihydrodiol epoxides. 



The importance of CYPIBI in chemical carcinogens is well 
illustrated in animal models in which metaboKtes of GYP 1 B 1 were 
shown to induce Prostate cancer risk [54,55]. GYP 1A2 is an 
important gene in catalyzing 2- and 4-hydroxylations of estrogens 
and metabolism of carcinogens. A major reason for the limited 
number of studies of heterocyclic amine (HGA) and cancer risk is 
the difficulty of assessing human exposure to HCAs. HCA 
concentrations depend on cooking methods and the "doneness" 
level of the meat or fish, hampering the development of a complete 
and standardized database of concentrations; any estimation of 
dietary intake from food-frequency questionnaires (FFQs) is thus 
likely to result in misclassification. Like other environmental 
chemical carcinogens, HCAs require metabolic activation by host 
enzymes to become genotoxic. Phase I enzymes, including 
cytochrome P450 1A2, can metabolically activate carcinogens to 
form genotoxic electrophUic intermediates [56]. The relative 
activity of these metabolizing enzymes, which is in large part 
genetically determined, is thought to be an important host 
determinant of cancer incidence. A number of epidemiological 
studies have evaluated the association between CYP1A2'''F, 
CYPIBI Leu432Val, Asn453Ser, and Arg48Gly polymorphisms 
and colorectal cancer risk, but the results remain inconclusive. In 
order to resolve this conflict, this meta-analysis of 39 eligible 
studies including 5,817 cases and 6,544 controls for GYP1A2*F 
(from 13 studies), 9,219 cases and 10,406 controls for CYPIBI 
Leu432Val (from 12 studies), 6,840 cases and 7,761 controls for 
CYPIBI Asn453Ser (from 8 studies), and 4,302 cases and 4,791 
controls for CYPIBI Arg48Gly (from 6 studies) was performed to 
derive a more precise estimation of the association between 
CYP1A2*F, CYPIBI Leu432Val, Asn453Ser, and Arg48Gly 
polymorphisms and risk of colorectal cancer. 

Overall, no significant association was found between CY- 
P1A2*F, CYPIBI Leu432Val, Asn453Ser, and Arg48Gly when 
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Figure 5. Begg's funnel plot of the meta-analysis of colorectal cancer risk and CYP1 A2*F polymorphism (homozygote model and 
dominant model). 

doi:1 0.1 371/journal.pone.01 00487.g005 



all the eligible studies were pooled into the meta-analysis. And in 
the subgroup, no evidence of significant association was also 
observed in any subgroup. Sachse et al. [33] in 2002 and Kiiry 
et al. [24] in 2007 reported that CYPIBI Leu432Val was not 
associated with increased the risk of colorectal cancer. Landi et al. 
[27] and Huber et al. [37] in 2005 reported that CYPIBI 
Leu432Val and Asn453Ser polymorphisms were also not associ- 
ated with increased the risk of colorectal cancer. Cleary et al. [18] 
in 2010 found that CYPIBI Leu432Val, Asn453Ser, and 
Arg48Gly were not associated with increased the risk of colorectal 
cancer. Sachse et al. [21] in 2002, Yoshida et al. [22] in 2007, Kiss 
et al. [23] in 2007, and Cleary et al. [18] reported that CYP1A2*F, 
was not associated with increased the risk of colorectal cancer. The 
results of our meta-analysis supported the negative association 
between CYP1A2*F, CYPIBI Leu432Val, Asn453Ser, and 
Arg48Gly polymorphisms and colorectal cancer risk. However, a 
careful matching should be considered in future larger genetic 
association studies including multiple ethnic groups. 

We noticed that 3 previous meta-analysis [33,57,58] had been 
reported on the colorectal cancer risk with CYP1A2*F, CYPIBI 
Leu432Val, and Asn453Ser polymorphisms. We have read with 



great interest the meta-analysis by Mei et al. [57] and Xie et al. 
[58]. Mei et al. [35] had 7 studies including 6,375 cases and 7,003 
controls. The pooled analysis suggested that no significant 
association was found between the CYPIBI Asn453Ser polymor- 
phism and the risk of colorectal cancer among Caucasians. Xie 
et al. [58] had 10 studies including 8,466 cases and 9,301 for 
Leu432Val. Their meta-analyses suggested that CYPIBI Leu432- 
Val were not associated with colorectal cancer risk. However, the 
study of Northwood et al. [30] should be excluded in the meta- 
analyses of Mei et al. [57] and Xie et al. [58] because they 
performed CYPIBI Leu432Val with colorectal adenoma risk but 
not colorectal cancer. Adopting the same search strategy as Mei 
et al. [57] and Xie et al. [58], we identified 4 additional eligible 
studies, which have not been included in the meta-analysis of Xie 
et al. [36]. Worthy of note, these 4 studies included 3,638 samples. 
Zhao et al. [33] included 1 1 studies. Their meta-analysis suggests 
that the CYP1A2*F polymorphism is a protective factor against 
CRC among Asians. The OR (95% CI) reported by Zhao et al. 
[33] for the study by Bae et al. [25] do not seem in line with the 
OR (95% CI) provided by Bae et al. [25] in their original 
publication. The OR (95% CI) reported by Zhao et al. [33] in 
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Figure 6. Begg's funnel plot of the meta-analysis of colorectal cancer risk and CYP1 B1 Leu432Val polymorphism (homozygote 
model and dominant model). 

doi:10.1371/journal.pone.0100487.g006 
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Figure 7. The Duval and Tweedie nonparametric "trim and fill" method's funnel plot funnel plot of the meta-analysis of colorectal 
cancer risk and CYP1B1 Arg48Gly polymorphism (additive model and dominant model). 

doi:1 0.1 371 /journal.pone.01 00487.g007 



additive model are 0.56 (0.38-0.84). Interestingly enough, after 
carefully studying the OR (95% CI) presented by Bae et al. [25], 
The OR (95% CI) were 1.77 (1.18-2.66). In addition, the study of 
Wang et al. [59] should be excluded in the meta-analysis of Zhao 
et al. [33] because the data on CYP1A2*F polymorphism with 
colorectal cancer risk did not be found in the study of Wang et al. 
[59]. Adopting the same search strategy as Zhao et al. [33], we 
identified 3 additional eligible studies, which have not been 
included in the meta-analysis of Zhao et al. [33]. Worthy of note, 
these 3 studies included 2687 samples. Having analyzed an almost 
twofold larger number of studies than the previous meta-analysis 
[33,57,58], our results seem to confu-m and establish the trend in 
the meta-analysis of CYP1A2*F, CYPIBI Leu432Val, Asn453Ser, 
and Arg48Gly polymorphisms that the data by the previous meta- 
analysis [33,57,58] had indicated. The results of the present meta- 
analysis are not in accordance with those reported by Zhao et al. 
[33]. Our meta-analysis indicates that CYP1A2*F are not 
associated with colorectal cancer risk. 

There are several limitations in this meta-analysis. First, the 
controls were not uniformly defined. Although most of them were 
common populations, some controls were population-based; other 



controls were hospital-based. Hence, non-differential misclassifi- 
cation bias is possible. Second, in the subgroup analysis may have 
had insufficient statistical power to check an association. Third, we 
were also unable to examine the interactions among gene- 
environment, lacking of the original data of the included studies 
limited our further evaluation of potential interactions, which may 
be an important component of the association between CY- 
P1A2*F, CYPIBI Leu432Val, Asn453Ser, and Arg48Gly poly- 
morphisms and environment and colorectal cancer risk. Last, our 
results were based on unadjusted published estimates. Because of 
data limitations, we were unable to adjust them such as age and 
alcohol consumption et al. 

In summary, this meta-analysis indicates that CYP1A2*F, 
CYPIBI Leu432Val, Asn453Ser, and Arg48Gly are not associ- 
ated with colorectal cancer. However, it is necessary to conduct 
large sample studies using standardized unbiased genotyping 
methods, homogeneous cancer patients and well-matched con- 
trols. Moreover, further studies estimating the effect of gene-gene 
and gene-environment interactions may eventually lead to our 
better, comprehensive understanding of the association between 
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Figure 8. The Duval and Tweedie nonparametric "trim and fill" method's funnel plot funnel plot of the meta-analysis of colorectal 
cancer risk and CYP1B1 Asn453Ser polymorphism (additive model and dominant model). 

doi:1 0.1 371/journal.pone.01 00487.g008 
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the CYP1A2*F, CYPIBI Leu432Val, Asn453Ser, and Arg48Gly 
polymorphisms and colorectal cancer risk. 
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